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Transition in Hypersonic Boundary Layers: Role of Dilatational Waves 
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Transition and turbulence production in a hypersonic boundary layer is investigated in a Mach 
6 quiet wind tunnel using Rayleigh-scattering visualization, fast-response pressure measurements, 
and particle image velocimetry. It is found that the second instability acoustic mode is the key 
modulator of the transition process. The second mode experiences a rapid growth and a very fast 
annihilation due to the effect of bulk viscosity. The second mode interacts strongly with the first 
vorticity mode to directly promote a fast growth of the latter and leads to immediate transition to 
turbulence. 
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Transition and turbulence production in hypersonic 
boundary layers have recently received considerable at¬ 
tention owing to their fundamental importance and 
strong relevance to the safety of hypersonic vehicle flight, 
including significant increases in aerodynamic heating, 
entropy production, and drag BE]. Compared to in¬ 
compressible flows 00, hypersonic transition is less un¬ 
derstood due to additional complexities, such as second 
and higher instability modes, and nonlinear coupling of 
different processes. The problem presents tremendous 
challenges to theory, experiment, and computation m\- 

During laminar-to-turbulence transition in compress¬ 
ible boundary layers, a second instability mode ap¬ 
pears at very high frequency and becomes stronger 
with increasing Mach number. Unlike the conventional 
first-mode instability associated with transverse vortical 
waves, the second-mode instabilities are acoustic, longi¬ 
tudinal waves reflected between the relative sonic line 
and the solid wall, and are of crucial importance at high 
Mach numbers 0i- Detailed work comparing the lin¬ 
ear growth of the second-mode with linear stability the¬ 
ory (LST) has been undertaken both in quiet and 

noisy tunnels [IMS]. 

Unfortunately, there exists conflicting understanding 
of the effects of the second-mode instability. Stetson m 
and Pruett et al. m believed that the second mode in¬ 
stability plays a dominant role in transition, while Li 
m, Bountin [18], and Dong et al. [19] took opposite 
view. Evidently, the dynamic process of the evolutions 
of the two modes and their interactions need to be clar¬ 
ified. Emanuel [20] and Gad-el-Hak [21], among others, 
recognized the important role of the second coefficient of 
viscosity in high-Mach-number flows [22] . By including a 
non-zero bulk viscosity as a correction term in a general, 
non-similar formulation of hypersonic laminar boundary 
layers, Emanuel m found a heat transfer rate well in 


excess of that based on Stokes’s hypothesis m predic¬ 
tion, i.e. assuming zero bulk viscosity. Despite these pro¬ 
gresses, neither qualitative visualization nor quantitative 
measurement of the complete transition process has yet 
been reported, and an underlying mechanism of tran¬ 
sition and its association to the bulk viscosity are still 
missing. 

In this Letter, we report an experimental study of 
the second-mode instability and its relevance to tur¬ 
bulence production. We use combined visualization- 
measurement techniques to obtain the most complete and 
detailed data set. We are able to display quantitatively 
the entire streamwise evolution of the first and second 
modes, from their appearance and growth to the transi¬ 
tion to turbulence. Through analysis, we reveal, for the 
first time, the underlying mechanism of the transition in 
hypersonic boundary layers and its strong connection to 
the second mode and the bulk viscosity. 

Our experimental model is a sharp-nosed flared cone, 
shown in Eig. Owing to the low freestream distur¬ 
bance in the quiet flow, natural transition does not occur 
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FIG. 1. Schematic of the model. 
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FIG. 2. Visualization of transition stages. Arrow A: appearance of second mode; B: second mode’s decay to almost zero; and 
C: onset of turbulence. The small vertical arrows indicate where a splice occurs. 


up to the end of an equivalent non-flared cone. There¬ 
fore, to obtain a whole picture of the transition process, a 
flared, instability-enhancing cone was studied under nat¬ 
ural freestream disturbance environment. The freestream 
Mach number 6. The angle of attack is zero. The mea¬ 
surements were carried out in a newly established Mach 6 
quiet wind tunnel (M6QT) at Peking University, which, 
at present, is one of three operational hypersonic quiet 
wind tunnels in the world m- A detailed description 
of the M6QT has been given by Zhang et al. [13]. The 
sharp-nosed flared cone has a straight 5° half-angle for 
the first 0.381/ of axial length, with a tangent flared re¬ 
gion of radius 3.581/ for the remaining 0.621/, designed 
to enhance the instability and thus favor transition. The 
total length of the cone is L = 260 mm. The origin of the 
coordinate system is located at the cone’s nose, with x 
being the streamwise coordinate along the cone’s surface, 
y being the coordinate normal to the cone’s surface, and 
2 : the transverse coordinate normal to the x-y plane. The 
freestream stagnation temperature and pressure are, re¬ 
spectively, 430 K and 0.9 MPa. The freestream velocity 
and unit Reynolds number are, respectively, 870 m/s and 
9.0 X 10^ per meter. 

In the first measurement reported here the whole tran¬ 
sition process is investigated using Rayleigh-scattering 
techniques, which can offer a clear but qualitative infor¬ 
mation on the evolution of the structures in the boundary 
layer [24] . In the Rayleigh-scattering images (RSI) shown 
in Fig.[^ independent measurements over four successive 
sections are spliced together to exhibit the whole evolu¬ 
tion of the second mode. The flow is from left to right. 
In the upstream region, the boundary layer is laminar. 
Different from incompressible flow [4], a linear instabil¬ 
ity of the hypersonic flow results in growing second-mode 
patterns. The unstable waves are weak at the early stage 
of the destabilized laminar flow, but eventually the am¬ 
plitude of these waves increases dramatically. The insta¬ 
bility waves are formed and persist for a long distance as 
shown in the region between arrow A and arrow B. The 
regular rope-like structures are second-mode waves and 
their wavelength is about 2 mm. After further nonlin¬ 
ear development, the second mode decays to small value 
as demonstrated in the region between arrows B and 
C, which will be quantitatively demonstrated in Fig. [^ 
The picture to the right of arrow C shows that a tur¬ 
bulent boundary layer appears immediately after transi¬ 


tion. The transition process is similar for different unit 
Reynolds numbers. Because the second-mode waves ex¬ 
ist for a long distance, a clear and complete visualization 
of the transition process has not been reported before. 

Surface-mounted PCB 132A31 fast-response piezoelec¬ 
tric pressure transducers were used to measure the pres¬ 
sure fluctuations on the model’s surface. The first sensor 
is located at x/L = 0.38, and the first four are positioned 
20 mm apart. From xjL = 0.62, the spacing between the 
neighboring sensors is 10 mm to acquire a detailed evolu¬ 
tion of the unstable disturbance. By analyzing the power 
spectral densities, we obtained detailed developments of 
frequencies and amplitudes of the first- and second-mode 
disturbances, as partially shown in Fig. Here the 
frequency results for the x-station locations down the 
model are shown with two distinct peaks, at 30 kHz and 
350 kHz. The first and second modes are clearly shown. 
At a frequency of 30 kHz, the first mode grows continu¬ 
ally to the onset of turbulence, while the second mode, 
at 350 kHz, grows initially and then decays very quickly, 
eventually reaching the relatively quiet region. The quick 



FIG. 3. Frequency spectrum at different x-positions. Both 
first mode and second mode are shown. 
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FIG. 4. (color online). PIV results of boundary layer development, (a) Velocity magnitude normalized with freestream velocity; 
(b) spanwise component of vorticity L)Jz\ (c) dilatation'd = dxU + dyV. Arrow A: appearance of second mode; B: second mode’s 
decay to almost zero; and C: onset of turbulence. 



FIG. 5. (color online). Viscosity dissipations induced by (a) vorticity; and (b) dilatation. Arrow A: appearance of second 
mode; B: second mode’s decay to almost zero; and G: onset of turbulence. 


growth of the second mode is consistent with Mack’s LST 
calculations [25] and experimental results obtained using 
hot-wire anemometry m- 

A second quantitative investigation of the transition 
process was carried out using a particle image velocime- 
try (PIV) system and PCB measurements simultane¬ 
ously. The RSI and PCB results show that injection of 
tracer particles in PIV may cause the transition to occur 
earlier, but the whole transition process is almost the 
same as shown in Fig. Based on this, a comparison 
between PCB data and PIV results can also be made. 
To solve the problem of very large in-plane displacement 
on PIV images, commonly encountered in hypersonic 
flows, we developed an improved image-preprocessing 
method m- Our strategy can extract the second-mode 
waves’ fine flow structures from PIV raw images, while 
traditional PIV methods are incapable of doing that. 
PIV images were at first rotated horizontally and pre- 
processed by the new method. The interrogation started 
with 256x256 pixels coarse samples, then 128x128 pix¬ 
els medium samples to capture the mean displacement 
fields from the time series. With the mean results as 
reference, instantaneous recordings were then evaluated 
with 64x64 pixels samples to resolve the fine structures 
in the boundary layer. Smaller samples are unavailable 
because of the low particle density near the wall. Instan¬ 
taneous PIV results of the boundary layer development 
are shown in Fig. where independent measurements 
over three successive sections are spliced together. 

We plot in Fig. the velocity magnitude normalized 
by the freestream velocity, the ^-component of vorticity. 


— dyU^ and dilatation = dxU + dyV. Note 
that if one writes u = U{y) u' with U{y) being the 
mean shear velocity, then there is ujz = —dyU{y) + 
where is the disturbance vorticity wave that is not 
shown in Fig. [^. Combined with Figs. and these 
figures provide a clear evolution process of the first and 
second modes. Notice a relatively quiet zone between 
arrows B and C with x/L G (0.60,0.62) in Figs. and 
[^ where the boundary layer thickness becomes thin with 
the decay of the second mode. At the same time, the 
viscous dissipation is at a low level (see Fig. [^. 

Of the new phenomena observed in our experiments, 
the most remarkable ones are the rapid annihilation of 
the second (dilatation) mode after its fast growth, and 
the associated very strong enhancement of the first (vor¬ 
tical) mode that leads to a quick transition. Although 
a few experiments have shown the dominance of the 
second-mode waves in the production of hypersonic tur¬ 
bulence and their disappearance in turbulent boundary 
layers m, the detailed evolution process of the second 
mode, especially its fast decay, has remained unclear. 
Schlieren or shadowgraph techniques can be used to ob¬ 
serve the second mode [T3|, but the fine flow structures 
in the boundary layer cannot be fully identified. There¬ 
fore, the above observed peculiar phenomena demand a 
physical interpretation, which we provide below. 

In hypersonic boundary layers the second mode domi¬ 
nates the transitional flow, see Fig. When the dilata¬ 
tion wave becomes sufficiently strong, it contributes con¬ 
siderably to the irreversible energy dissipation (Figj^. 
Indeed, it is well-known that, in addition to shear viscos- 
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ity the compressibility of a viscous flow with nonzero 
dilatation brings in a longitudinal viscosity /r' = + |/i, 

where is the bulk viscosity, so that the pressure p is 
modified to p — Accordingly, the dissipation func¬ 
tion 4> per unit volume can be explicitly split to the con¬ 
tributions of vorticity and dilatation [221 EH, where a 
divergence term of undetermined sign is neglected: 

4 > =; uj = \lv\, ( 1 ) 

where p is the shear viscosity. In monatomic gases, 
Stokes’s hypothesis /i^ = 0 holds, implying p' = |/i. 
But in diatomic and polyatomic gases /i^ 7 ^ 0. Both p 
and fii) are temperature dependent. Figure shows the 
shear-induced dissipation which is largely indepen¬ 
dent of the compressibility but affected by temperature. 
The term shown in Fig. is a pure compressibil¬ 
ity effect and represents the dilatation-induced dissipa¬ 
tion. Here, a remarkable fact is that the bulk viscosity 
is not only a parameter of material property but also a 
dispersive function depending on frequency. While in a 
non-moving gas = 0 (/i), it can be magnified signifi¬ 
cantly at high-frequency circumstances [28]. The afore¬ 
mentioned high-frequency trapped acoustic waves rele¬ 
vant to hypersonic second-mode instability just provide 
such a circumstance. In our study, was chosen as 
0.73/i [ 29 ], its value in quiet ^^2 at T = 293 K. Given the 
high frequency of 350 kHz and higher temperature, the 
value we use is certainly an over-conservative estimate 
and very likely to be much larger in hypersonic bound¬ 
ary layers (could be by a factor of 10^ higher). 

Now, the amplitude of the instability waves and the 
viscous dissipation can be calculated by combining the 
simultaneous pressure and PIV measurements. Figure 
displays more clearly the peak values of /icj^ and 
as well as their correlation with the spatial evolutions 
of first- and second-instability-wave amplitudes. Fig¬ 
ure [Slshows the viscosity dissipation induced by vorticity 
(Fig7[^) and dissipation (Fig. [^). The viscous dissipa¬ 
tion values were extracted from Figs. and [^. The 
dissipation caused by was sampled at the local peak 
points of the region where 'd 7 ^ 0. Each quantity was 
normalized by its maximum value. 

Physically, the aforementioned modal interactions can 
be explained by the strong (u;,'d)-couplings that exist 
generically both in the interior of the boundary layer and 
at the wall [26| . The internal coupling occurs via the com¬ 
pressible vorticity transport equation. In the flow region 
away from the viscous sublayer with negligible normal 
motion, the disturbance vorticity is governed by 


duj' .duj' ^ , 


( 2 ) 


Our experiment has found that in the range of x/L G 
(0.47,0.53), the vortical and acoustic waves have almost 



FIG. 6. (color online). Comparison of wave amplitude and 
dissipation function by shear and dilatation. Amplitude of 
first mode is taken from PCB, while that of second mode from 
both PCB and PIV. Dissipation caused by dilatation is cal¬ 
culated for both pb = 0.73/i (an over-conservative estimate) 
and ph = 0. 


the same phase speed c ISO], so both can be approxi¬ 
mately expressed as functions of 7 / = x — ct, with c 7 ^ ii 
in this region. Thus, (§ yields the following coupling 
relation between the phase-locked traveling Uz-wsMe and 
'd-wave 


( 3 ) 

indicating that the cj-wave can be significantly magni¬ 
fied by the 'd-wave at its peak-value zone and near the 
critical layer. The theory of the phase-locked nonlinear 
interaction has been well developed, known as nonlin¬ 
ear critical-layer analysis (see the review by (SB). It 
states that the second mode will significantly prompt the 
growth of the first mode. Our experimental results show 
that the small-amplitude first mode has gone through 
a remarkable growth while the large-amplitude second 
mode decays sharply, indicating that the energy has been 
transferred from the high-frequency disturbances to the 
low-frequency disturbances. This can be viewed as a di¬ 
rect validation of the theory. 

On the other hand, the boundary (co’,p)-coupling oc¬ 
curs via viscosity and the no-slip condition, which can be 
estimated more accurately than Eq. (3). Consider the in¬ 
teraction of two traveling waves of the form off] = x — ct: 
the normal-stress wave H = p — = Ii(kQr]) and the 

shear-stress wave with ke 1 and k^j kg be¬ 

ing the wave numbers of second and first modes, respec- 
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tively. The vorticity creation rate at the wall is measured 
by cr = —vdyUJ dX y = 0[34]. Then due to the no-slip 
condition, the on-wall viscous (co’,p)-couphng takes the 
form[29,34] 

1 f)TJ 

a{x,t) =-ly—= — = kell'iri) at y = 0, (4) 

indicating that the high-frequency 'd-wave must create 
new cj-wave of the same resulting in a high-frequency 
Stokes layer, which directly prepares the appearance of 
strong turbulent fluctuations. The extra large bulk vis¬ 
cosity /i 5 must significantly enhance this boundary cou¬ 
pling. 

With these (c(;,p)-couplings inside the flow and at the 
wall, why the first-mode instability with small vorticity 
wave number kg is quickly enhanced and gains 

high-frequency components right after the dissipation of 
second mode can now be understood. It may be men¬ 
tioned that symmetric to Eq. (1) there is 

^ = vk^uj'{-q) at y = 0, (5) 

which represents how the low-frequency first-mode wave 
alters the Neumann condition of the dilatation equation. 
But since vk^^ ^ Re~^^‘^k^ ^ 1, this effect is pretty weak 
or even negligible. This confirms the preceding analysis 
on the one-way interaction of fast-mode and slow-mode 
instability waves. 

In summary, our results suggest the following route 
from second-mode instability to turbulence production. 
At first, the free-stream disturbances penetrate into the 
boundary layer near the leading edge of the cone and 
then evolves as the linear stability theory predicts. The 
unstable waves (the first mode and the second mode) 
will be amplified. When the amplitudes of these unsta¬ 
ble become large enough, the nonlinear interactions may 
be triggered. The nonlinear coupling of dilatation waves 
and vortical waves inside the boundary layer and their 
linear coupling at the wall enable the highly localized 
strong peak of dilatation wave to enhance the first-mode 
instability and thereby cause the production of turbu¬ 
lence. In particular, the bulk viscosity /i^, which has so 
far been overlooked by most theoretical studies and nu¬ 
merical simulations, serves as a strong magnifier not only 
in the very strong dissipation peak of high-frequency di¬ 
latation waves but also in the viscous (u?,'d)-coupling at 
the wall. The proposed scenario may inspire new strate¬ 
gies for the control of transition in hypersonic boundary 
layers. 
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